A Minimum Principle in 
Codon-Anticodon Interaction 

A. Sciarrino''''^'*, P. Sorba^ 

" Dipartimento di Scienze Fisiche, Universita di Napoli ''Federico //" 
■ Complesso Universitario di Monte S. Angela 

CN ; Via Cinthia, 1-80126 Napoli, Italy 

''I.N.F.N., Sezione di Napol 
sciarrino@na. infn. it 



o 
O 



Laboratoire d'Annecy-le-Vieux de Physique Theorique LAPTH 
CNRS, UMR 5108, associee a I'Universite de Savoie 
BP 110, F-74941 Annecy-le-Vieux Cedex, France 
sorba@lapp.in2p3.fr 

6 

X) 

D^i Abstract 

Imposing a minimum principle in the framework of the so called crystal basis model of the genetic 
>- ' code, we determine the structure of the minimum set of anticodons which allows the translational- 
^ ■ transcription for animal mitochondrial code. The results are in very good agreement with the ob- 
! served anticodons. 
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1. Introduction 



The translational-transcription process from DNA to proteins is a very complex process carried 
on in several steps. A key step is the translation from coding sequences of nucleotides in mRNA to 
the proteins chaines. In this process a role is played by the tRNA in which a triplet of nucleotides 
(anticodon) pairs to the triplet of nucleotides (codon) reading the genetic information. Since there 
are 60 codons (in mitochondrial code) specifying amino acids, the cell should contain 60 different 
tRNA molecules, each with a different anticodon in order to have a pairing codon anticodon following 
the usual Watson-Crick pattern, i.d. the pairing respectively between the nucleotides C and G, and 
U and A. Actually, however, the number of observed anticodons is less than 60. This implies that an 
anticodon may pair to more than one codon. Already in the middle of the sixties , it wa s realized that 
the pairing anticodon- codon does not follow the standard rule and Crick (j Crick I . Il966l ) proposed, on 
the basis of the base-pair stereochemistry, the "wobble hypothesis". According to this hypothesis a 
single tRNA type, with a a specified anticodon, is able to recognize two or more codons in particular 
differing only in the third nucleotide, i.e only the first two nucleotides of a codon triplet in mRNA 
have the standard precise pairing with the bases of the tRNA anticodon while the first nucleotide in 
the anticodon may pair to more than a nucleotide in the third position of the codon. 

This rule has been subsequently widely confirmed and extended, Av ith a better understanding 
of the chemical nucleotide modifi c ation s, for a review see (lAgris 1 . 1200J). Since the years seventies 
the questions were raised (jjukes 1 . 119771 ): how many anticodons do we need? which anticodons do 
manifest? 

In order to explain which anticodon do manifest two main hypothesis have been advanced: 

1. The conventional wobble versatility hypothesis assumes that the the first position of anticodon 
should have G (U) to read for codon with Y (respectively R) in third position. 

2. The codon adaptation hypothesis states that the first position of anticodon should pair the 
most abondant codon in the family of synonymous codons. 

For a comparison and discussion of the t wo hypothesis in fungal mito chondrial genomes and fo r 
marine bivalve mitochondrial genomes, see (jCaruUo and Xia I . l2000l ) and (jHong Yu and Qi li 1 . 120 111 ). 

In order to have a correct translation process between codons and amino-acids in the mitochon- 
drial code we need a minimum number of 22 anticodons. In fact, in this code, the 20 amino-acids 
(a. a) are encoded by 2 sextets, 6 quadruplets and 12 doublets of codons. Considering a sextet as the 
sum of a quadruplet and a doublet, we need to dispose at least of 22 anti-codons, of which 8 should 
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200.1 IWilhelm and Nikola jewa I . l2004l : iNikolajewa et al. 



The data seem to confirm the empirical rule that the most 
used anticodons have as second and third nucleotide, respectively, the complementary to the first 
and second nucleotide of the codons, while the first nucleotide is U for the anticodons pairing the 
quadruplets, G and U for the anticodons pairing, respectively, the doublets ending with a pyrimidine 
and with a purine, with exception of Met. 
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The aim of this paper is to propos e a mathematical ap proach in the framework of the "crystal 
basis model" model of the genetic code (IFrappat et all . ll998l ). to determine which anticodon is chosen 
to translate the genetic information stored into the quadruplets and the doublets of codons. More 
generally, the idea is to require the minimization of a suitable operator or function, mathematically 
expressed in terms of the quantities defined in the model, to explain why and which anticodon is 
used to "read" more than a codoE0. 



Let us very quickly recall the main ideas of the model introduced in (IFrappat et al I . Il998l ). 



for a review and some applications see flFrappat et al I . l200ll ). In that paper we have proposed a 



mathematical framework in which the codons appear as composite states of nucleotides. The four 
nucleotides being assigned to the fundamental irreducible representation (irrep.) of the quantum 
group lAq{su{2) ®su{2)) in the limit g — > 0, the codons are obtained as tensor product of nucleotides. 
Indeed, the properties of quantum group representations in the limit g — > 0, or crystal basis, are 
crucial to take into account the fact that a codon is an ordered triple of nucleotides. The nucleotide 
content of the (|, |) (fundamental) representation of Uq^o{su{2) © su{2)), i.e. the eigenvalues of 
Jh3, Jv3y is chosen as follows: 
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+ : 



U 



1 +h 



G 
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A 



1. 
'2' 



(1) 



where the first su{2) - denoted su{2)h- corresponds to the distinction between the purine bases A,G 
and the pyrimidine ones C,U and the second one - denoted su{2)v - corresponds to the complemen- 
tarity rule C/G and U/A, Thus to represent a codon, we have to perform the tensor product of three 
|) or fundamental representations oilAq^Q{su{2)®su{2)) and we get the results, reported in Table 



'2' 2 



El where we have als o written the observed anticodon for the mitochondria of animals taken from 



( ISprinzl et al. I . Il998[ ). Really in the present paper we use the minimum principle in a reduced form 



as we are only interested to find the composition of the minimum number of anticodons. However 
in the last section we hint at some more general application of our schema. 

The paper is organised as follows: in Sec. 2 we present the minimum principle, in Sec. 3 we apply 
the prin ciple to the mitochon drial code for animals and we compare our theoretical results with the 



data of f lSprinzl et al. I . Il998l ). In the final Section we give a short summary and some highlights on 



future developments and applications. 

2. The "minimum" principle 

Given a codoni XY Z {X, Y, Z e {C, A, G, U}) we conjecture that an anticodon X^-Y^Z", where 
ya^a _ YcXc, Nc denoting the nucleotide complementary to the nucleotide N according to the 
Watson-Crick pairing rul^, pairs to the codon XYZ, i.e. it is most used to "read" the codon XYZ 



We do not discuss here the chemical modified structm-c of the nucleotides, e.g see (lAgrisl . l2004l) . 
^In the paper we use the notation N = C, A,G,U.; R ^ G, A. (purine); Y = C,U. (pyrimidine). 
•^This property is observed to be verified in most, but not in all, the observed cases. To simplify we shall assume it. 
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if it minimizes the operator T, explicitly written in eq. ([2]) and computed between the "states" , which 
can be read from Table Ej describing the codon and anticodon in the "crystal basis model". We write 
both codons (c) and anticodons (a) in 5" — > 3" direction. As an anticodon is antiparallel to codon, 
the 1st nucleotide (respectively the 3rd nucleotide) of the anticodon is paired to the 3rd (respectively 
the 1st) nucleotide of the codon. 

r = 8cH 4 ■ 4 + 8cy 4 ■ 4 (2) 

where: 

• ch-Cv are constants depending on the "biological species" and weakly depending on the encoded 
a.a., as we will later specify. 



Jjilf, J y (resp. J^, Jy) are the labels of Ug^o{su{2) H®su{2)y) specifying the state flFrappat et al I . 



19981 ) describing the codon XYZ (resp. the anticodon NYcXc pairing the codon XYZ). 



J^- 4i (a = H, V) should be read as 

4-4 = l[{jl® Jc^y - (f^y - (^1)'} (3) 

and Ja® Ja = Jql stands for the irreducible representation which the codon-anticodon state 
under consi deration belongs t o, the tensor product of and being performed according to 



the rule of f lKashiwara 1 . 11990[ ). choosing the codon as first vector and the anticodon as second 

^ 2 

vector. Note that Ja should be read as the Casimir operator whose eigenvalues are given by 

Ja{Ja + 1). 

For example the value of T between the anticodon UUU and the codon AAC is, using Table H) 

<UUU\T\AAC >= -QcH + lScv (4) 

As we are interested in finding the composition of the 22 anticodons, minimun number to ensure 
a faithful translation, we shall assume that the used anticodon for each quartet and each doublet 
is the one which minimizes the averaged value of the operator given in eq.(|2]), the average being 
performed over the 4 (2) codons for quadruplets (doublets), see next section. Indeed it is well known 
that synonymous codons are not used with equal frequency. Therefore, in finding the structure of the 
anticodons in the minimum set, it appears reasonable that the codon usage probability plays a role 
in the determination of the chosen anticodon. If the codon XY Z is used more frequently than the 
codon XYZ' , the codon XY Z should give an impact larger than the codon XYZ' in determining the 
choice of the anticodon in the T. Therefore, Tav appears more appropriate than T for our purpose. 



3. Structure of the minimum number of anticodons 

According to our conjecture on the existence of a minimum principle we determine, for each 
quadruplet (q) and each doublet (d), the anticodon which minimizes the averaged value Tav of the 
operator T (see below). We analyse separately the case of quadruplets and doublets. 
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3.1. Quadruplets 

Let us give an example of what we mean by averaged value of T. For example let us consider the 
anticodon CAC for the a. a. Val, we have to compute 

Tav{CAC,Val) = ^ < CAC\r\GUN > 

N 

= P^c < CAC\T\GUC > < CAC\T\GUU > 
< CAC\T\GUG > +Pl < GAG\r\GUA > 
= 2{P^ + p^ + p^ + PI) CH + (6P^ + + + 2Pl) cv 
= 2cH + [QP^ + 2{l-P^)]cv (5) 

In the computation we have to take into account the codon usage frequency or relative percentage 
of the appearance of each codon in the quadruplet and we have denoted with the codon usage 
frequency for codon ending with N. Really we need to introduce the following four positive frequencies 
Py P^, P|,P^, with the normalization condition: 

P^ + P'r = PI + P^ = 1 (6) 

where, respectively, Py, P^, P| and P^ denote the relative usage frequency of the codons ending 
with nucleotides C,U (pyrimidine) , G,A (purine), C,G and U,A. From Table [3] we can compute the 
value which we report in Table El 

3.2. Doublets 

In the computation we have to take into account the codon usage frequency in the doublet. Now 
we need to introduce the following four positive frequencies P^, P^, P^, P^, with the normalization 
condition 

P^ + P^ = Pj + P^ = 1 (7) 

As example let us compute the averaged value of T for Asp. Let us consider the anticodon CUC we 
have to compute 

rav{CUG,Asp) = ^ P^ < CUG\r\GAY >= P^ < GUC\T\GAG > +P^ < CUG\r\GAU > 

Y 

= 2 ch + 18 Cv (8) 

From Table HI we can compute the values, which we report in Table El 
Let us remark that: 

• for all a.a. the contribution of suy{2) verifies the same property than for the quadruplets and, 
moreover, is not depending on the codon usage; 

• for 4 a.a. the contribution of suy[2) is the same for all anticodon. 
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From the above remarks we easily realize that the case of doublets is more complicated than the 
one of the quadruplets. In some sense the contribution of sny(2) plays a role only in establishing 
the most preferred anticodon. Moreover, as we do not want a priori to exclude any anticodon, we 
have to face the possibility that an anticodon can be chosen to read for more than one doublet. In 
order to avoid this problem, in contradiction with the requirement of a faithful translation process, 
we make the following choice: 

1. the sign of the constant ch for the doublets ending with a purine is the opposite of the sign of 
the doublets ending with a pyrimidine with the same dinucleotide (if it does exist )0. 

2. the sign of ch for the 8 weak dinucleotides encoding doublets is positive for the following 4 
doublets UUY, UAY, AUY, AAY and negative for the remaining 4, i.e. CAY, UGY, AGY, 
GAY. 

and fix the following procedure, while considering doublets with the same dinucleotide: 

1. first we select, among the four possible anticodons, the one giving the lowest value for T 
averaged on the two codons of each doublet and assign this anticodon to the corresponding 
doublet. 

2. then the anticodon reading the second doublet is chosen between the two ones containing as 
a first nucleotide a purine, resp. a pyrimidine, if the first nucleotide of the anticodon already 
determined for the first doublet is a pyrimidine (resp. a purine). 

As an illustration, we take the case of the Cys and Trp amino acids. The anticodon GCA can 
minimize both of them, but more Cys (due to the value —6 for the cy coefficient) than Trp (with 
value 2 for the same cy coefficient). Thus GCA will be taken as the anticodon relative to Cys, while 
the choice for the Trp anticodon will be made between UCA and CCA, that is the two candidates 
with a pyrimidine as a first nucleotide, the anticodon GCA starting with a purine. 

Let us remark that, even if the above assumptions seem rather ad hoc, indeed a general symmetric 
pattern shows up: for half of a. a. ch is positive and for the other half is negative; the first set of 4 
dinucleotides involves only 'weak" nucleotides, the second one a "strong" nucleotide; the dinucleotides 
XY and YX correspond to the same sign. 

3.3. Discussion 

For all quadruplets, from Tabled we remark that for ch > and cy < the anticodon minimizing 
the average value of T has the composition UXcYc. For Leu, Val and Thr a very weak condition for 
codon usage frequency has to be satisfied, i.e for the first two a.a. F| > 0,25 and for the last one 
Py > 0, 12£@. The results are in agreement with the observed anticodons, see ( Sprinzl et al. . 1998 ) 
and Table El 



■*Wc call dinucleotide the first two nucleotids of the codon. 

^ The constraint on the codon usage frequency can be released, imposing a suitable condition between ch and cy. 
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For doublets, we remark that, with the choice of sign of ch above specified and cy > for all a.a., 
the anticodons minimizing the average value of T are in agreement with the observed anticodon, see 
(jSprinzl et al. I . Il998l ) and Table [21 We summarize in Table [T] the results for the doublets. 
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Tabic 1: Anticodon minimizing the operator T, averaged over the two codons, for any amino acid encoded by a 
doublet, specifying the sign of ch- 



Let US remark that we find that for Met the anticodon is not UAU, as it should be expected frorn 



the e mpirical rule above quoted, but CAU which seems in agrement with the data, see (ISprinzl et al. 
1998h . 



4. Conclusions 

We have found that the anticodons minimizing the conjectured operator T given in eq.([2]), aver- 
aged over the concerned multiplets, are in very good agreement, the results depending only on the 
signs of the two coupling constants, with the observed ones, even if we have made comparison with 
a limited database. 

The fact that the crystal basis model is able to explain, in a relatively simple way, the ob- 
served anticodon-codon pairing which has its roots on the stereochemical properties of nucleotides 



( iLim and Curran I . 1200 ll ) strongly suggests that our modeling is able to incorporate some crucial 
features of the complex physico-chemical structure of the genetic code. 

It is rather clear that the operator T can be looked at as a codon-anticodon interaction operator. 
In this spirit, given a codon, the selected anti-codon for fixing the corresponding amino- acid, is 
determined as the one which minimizes the interaction. It is hard at this point to be more specific 
on the physico-chemical aspects of this quantity. 
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Indeed, the analysis in ( iPercudani and Ottonello I . 119991 ) suggests that the use of the wobble 
behaviour i s also dictated from the requirement of the optimization of translational efficiency and in 
( iLehmann 1 . 12000| ) the free-energy change of anticodon-codon interaction has been put into relation 
wit h the dissociation time, depending on the molecular structure of the encoded amino-acid, while 
in ( iLim and Curran I . l200ll ) emphasis has been put in the hydrogen and ionic bondings. 

It is intriguing that complex behaviour involving thermodynamical considerations as well as 
evolutionary effects can be cast in a single simple mininum principle. 

It might be interesting to note, among the different previously obtained applications of the crystal 
basis model, that this model has previously allowed to establish a pattern of correlations between 
the physico-chemical p roperties of the amino-acids and the assignment of the corresponding coding 
codons in the model ( iFrappat et al I . 120021 ) . Incidentally let us rernark th at the model explains 
the symmetry codon anticodon remarked in (jWilhelm and Nikolajewa I . I2004I ). Let us stress that our 
modeling has a very peculiar feature which makes it very different from the standard 4- letter alphabet, 
used to identify the nucleotides, as well as with the usual modeling of nucleotide chain as spin chain. 
Indeed the identification of the nucleotides with the fundamental irrep. of Uq{su{2)H © su{2)y) 
introduces a sort of double "bio-spin", which allows the description of any ordered sequence of n 
nucleotides as as state of an irrep. and allows to describe interactions using the standard powerful 
mathematical language used in physical spin models. 

In the present paper we have faced the problem to find the structure of the mimimum set of 
anticodons and, then, we have used a very simple form for the operator T, with the main aim to 
present a simple, mathematical modeling of the extremely complex codon-anticodon interaction. We 
have not at all discussed the possible appearance of any other anticodon, which should require a 
more quantitative discussion. For such analysis, as well as for the eukaryotic code, the situation may 
be different and more than an anticodon may pair to a quartet. For this future aim we have here 
reported Tables [3] and HI 

The pattern, which in the general case may show up, is undoubtedly more complicated, depending 
on the biological species and on the concerned biosynthesis process, but it is natural to argue that 
the usage of anticodons exhibits the general feature to assure an "efficient" translation process by 
a number of anticodons, minimum with respect to the involved constraints. A more refined and 
quantitative analysis, as well as comparison with other organisms, which should require more data, 
depends on the value of these constants. Very likely, it might happen that the assumption of the 
"universal" feature of ch and cy should be released and that the expression of the operator eq.([2]) 
should be modified, for example by adding a term of "spin-spin" interaction of the type 



^9h Jh,3'^h,3 + Jvs'^v 



V,3 



(9) 



where the values of Jh,3 and Jy^s, both for codons and anticodons, can be read out from their 
nucleotide composition, see Table |2j However the pattern which shows up in Tables |3] and HI with 
the values of the coefficients equal in pairs, strongly suggests that the minimum number of anticodons 
should be 32 (3 for the sextets, 2 for quadruplets and triplet and 1 for doublets and singlets). 
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Table 2: The vertebral mitochondrial code. The upper lafel denotes differen t irreducible representations. We list the 
most used anticodons for mitochondria of animals, see (jSprinzl et al. I . I1998I ). In bold- red (italic- blue) the anticodons 
reading quadruplets (resp. doublets). 



a. a 


CO don 


Kh,c 


Kv,c 




Kv,u 


Kh,g 




Kh,a 


Kv,A 


Pro 


CCC 


18 


-10 


-6 


-10 


18 


-30 


-6 


-30 




ecu 


6 


-10 


-10 


-10 


6 


-30 


-10 


-30 




CCG 


18 


-6 


-6 


-6 


18 


-10 


-6 


-10 




CCA 


6 


-6 


-10 


-6 


6 


-10 


-10 


-10 


Leu 


cue 


2 


-10 


-10 


-10 


2 


-30 


-10 


-30 




cuu 


2 


-10 


6 


-10 


2 


-30 


6 


-30 




cue 


2 


-6 


-10 


-6 


2 


-10 


-10 


-10 




CUA 


2 


-6 


6 


-6 


2 


-10 


6 


-10 


Arg 


CGC 


18 


2 


-6 


2 


18 


-6 


-6 


-6 




ecu 


6 


2 


-10 


2 


6 


-6 


-10 


-6 




CGG 


18 


2 


-6 


2 


18 


2 


-6 


2 




CGA 


6 


2 


-10 


2 


6 


2 


-10 


2 


Ala 


GCC 


18 


6 


-6 


6 


18 


-22 


-6 


-22 




ecu 


6 


6 


-10 


6 


6 


-22 


-10 


-22 




GCG 


18 


2 


-6 


2 


18 


6 


-6 


6 




GCA 


6 


2 


-10 


2 


6 


6 


-10 


6 


Gly 


GGC 


18 


18 


-6 


18 


18 


-6 


-6 


-6 




ecu 


6 


18 


-10 


18 


6 


-6 


-10 


-6 




CGG 


18 


18 


-6 


18 


18 


18 


-6 


18 




GGA 


6 


18 


-10 


18 


6 


18 


-10 


18 


Val 


cue 


2 


6 


-10 


6 


2 


-22 


-10 


-22 




GUU 


2 


6 


6 


6 


2 


-22 


6 


-22 




GUG 


2 


2 


-10 


2 


2 


6 


-10 


6 




GUA 


2 


2 


6 


2 


2 


6 


6 


6 


Ser 


UCC 


6 


-10 


-10 


-10 


6 


-30 


-10 


-30 




UCU 


2 


-10 


2 


-10 


2 


-30 


2 


-30 




UCG 


6 


-6 


-10 


-6 


6 


-10 


-10 


-10 




UCA 


2 


-6 


2 


-6 


2 


-10 


2 


-10 


Thr 


ACC 


6 


6 


-10 


6 


6 


-22 


-10 


-22 




ACU 


2 


6 


2 


6 


2 


-22 


2 


-22 




ACG 


6 


2 


-10 


2 


6 


6 


-10 


6 




ACA 


2 


2 


2 


2 


2 


6 


2 


6 



Tabic 3: Values of the coefficient multiplying ch [Kh = 8J|f • Jfj) and cy {Ky ~ 8Jy • Jy) computed from the value 
of the tensor product of the codon XY Z with the anticodon NYcXc, for the quadruplets. 
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a. a 


codon 




^ V ,o 


Khu 


Kvii 

^ V ,U 


Kmc 


^ V ,ljr 


Km A 


KvA 


His 


CAC 


2 


2 


-10 


2 


2 


-6 


-10 


-6 




CAT! 


2 


2 




2 


2 


-fi 




-fi 


Gin 


GAG 


2 


2 


-10 


2 


2 


2 


6 


2 




GA A 


2 


2 




2 


2 


2 




2 


Phe 


UUG 


-10 


-10 


-6 


-10 


-10 


-30 


-6 


-30 




UUU 


6 


-10 


18 


-10 


6 


-30 


18 


-30 


Leu 


UUG 


-10 


-6 


-6 


-6 


-10 


-10 


-6 


-10 




UUA 


6 


-6 


18 


-6 


6 


-10 


18 


-10 


Cys 


UGG 


6 


2 


-10 


2 


6 


-6 


-10 


-6 




UGU 


2 


2 


2 


2 


2 


-6 


2 


-6 


Trp 


UGG 


6 


2 


-10 


2 


6 


2 


-10 


2 




UGA 


2 


2 


2 


2 


2 


2 


2 


2 


Tyr 


UAG 


-10 


2 


-6 


2 


-10 


-6 


-6 


-6 




UAU 


6 


2 


18 


2 


6 


-6 


18 


-6 


Asp 


GAG 


2 


18 


-10 


18 


2 


-6 


-10 


-6 




GAU 


2 


18 


6 


18 


2 


-6 


6 


-6 


Glu 


GAG 


2 


18 


-10 


18 


2 


18 


-10 


18 




GAA 


2 


18 


6 


18 


2 


18 


6 


18 


He 


AUG 


-10 


6 


-6 


6 


-10 


-22 


-6 


-22 




AUU 


6 


6 


18 


6 


6 


-22 


18 


-22 


Met 


AUG 


-10 


2 


-6 


2 


-10 


6 


-6 


6 




AUA 


6 


2 


18 


2 


6 


6 


18 


6 


Ser 


AGG 


6 


18 


-10 


18 


6 


-6 


-10 


-6 




AGU 


2 


18 


2 


18 


2 


-6 


2 


-6 


Asn 


AAG 


-10 


18 


-6 


18 


-10 


-6 


-6 


-6 




AAU 


6 


18 


18 


18 


6 


-6 


18 


-6 


Lys 


AAG 


-10 


18 


-6 


18 


-10 


18 


-6 


18 




AAA 


6 


18 


18 


18 


6 


18 


18 


18 



Table 4: Values of the coefRcient multiplying ch [Kh = 8J|f • Jfj) and cv {Ky = 8Jy • Jy) computed from the value 
of the tensor product of the codon XY Z with the anticodon NYcXc, for the doublets. 
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a. a 


anticodon 


coeff. Ch 


coeff. Cy 


Pro 


CGG 
UGG 
GGG 
AGG 


18P| +6(1-P|) 
-6P| -10(1-P|) 
ISPJ +6(1-PJ) 
-6P| -10(1-P| 


-lOP^ - 6(1-P^) 
-lOP^ - 6(1-P^) 
-30P^ - lO(l-P^) 
-30P^ - lO(l-P^) 


Leu 


CAG 
UAG 
GAG 
AAG 


2 

-10P| + 6(1-P|) 
2 

-10P| + 6(1-PJ) 


-lOP^ - 6(1-P^) 
-lOP^ - 6(1-P^) 
-30P^ - lO(l-P^) 
-30P^ - lO(l-P^) 


Arg 


CGG 
UGG 
GGG 
AGG 


18PJ +6(1-P^) 
-6P| -10(1-P|) 
ISPJ +6(1-P|) 
-6P| - 10(1-P|) 


2 
2 

-6P^ +2(1-P^) 
-6P^ +2(1-P^) 


Ala 


GGG 
UGG 
GGG 
AGG 


18PJ + 6(1-PJ) 
-6P| - 10(1-P|) 
ISPJ + 6(1-P|) 
-6PJ - 10(1-P|) 


6P^ + 2(1-P^) 
6P^ + 2(1-P^) 
-22P^ + 6(1-P^) 
-22P^ + 6(1-P^) 


Gly 


GGG 
UGG 
GGG 
AGG 


18PJ + 6(1-PJ) 
-6PJ -10(1-P|) 
18P| + 6(1-P|) 
-6P| -10(1-P|) 


18 
18 

6P^ + 18(1-P^) 
6P-^ + 18(1-P^) 


Val 


GAG 
UAG 
GAG 
AAG 


2 

-lOPj + 6(1-P|) 
2 

-10P| + 6(1-P|) 


QP^ + 2(1-P^) 
6P^ + 18(1-P^) 
-22P^ + 6(1-P^) 
-22P^ + 6(1-P^) 


Ser 


GGA 
UGA 
GGA 
AGA 


6P| + 2(1-PJ) 
-10P| + 2(1-PJ) 

6P| + 2(1-P|) 
-10P| + 2(1-P|) 


-lOP^ - 6(1-P^) 
-lOP^ - 6(1-P^.) 
-30P^ - lO(l-P^) 

-30P^ - io(i-p;^) 


Thr 


GGU 
UGU 
GGU 
AGU 


6P| +2(1-P|) 
-10P| +2(1-PJ) 

6P| +2(1-P|) 
-10P| +2(1-P|) 


6P^ + 2(1-P^) 
6P^ + 2(1-P^) 
-22P^ + 6(1-P-^) 
-22P^ + 6(1-P^) 



Table 5: Value of the coefficients multiplying ch and cy in Tav , computed for any anticodon and averaged over the 
four codons for each quadruplet. 
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cl.cL 


anticodon 


coeff. ch 


coeff. Cy 


His 


CUG 
UUG 
GUG 
AUG 


2 

-IQP^ + 6(1-P^) 
2 

-IQP^ + 6(1-P^) 


2 
2 
-6 
-6 


Gin 


CUG 
UUG 
GUG 
AUG 


2 

-l^P^ + 6(1-P^) 
2 
6 


2 
2 
2 
2 


Phe 


CAA 
UAA 
GAA 
AAA 


-IQP^ + 6(1-P^) 
-QP^ + 18(1-P^) 
-10P($^ + 6(1-P^) 
-6P^ + 18(1-P^) 


-10 
-10 
-30 
-30 


Leu 


CAA 
UAA 
GAA 
AAA 


-lOP^ + 6(1-P^) 
-QP^ + ls\l-P^) 
-W^ + 6(1-P^) 
-6Pj + 18(l-Pj) 


-6 
-6 
-10 
-10 


Cys 


CCA 
UCA 
GCA 
ACA 


6P^ + 2(1-P^) 
-IQP^ + 2(1-P^) 

6P^ + 2(1-P^) 
-IQP^ + 2(1-P^) 


2 
2 
-6 
-6 


Trp 


CCA 
UCA 
GCA 
ACA 


6P^ + 2(1-P<^) 
-10P($^ + 2(1-Pj) 

6P^ + 2(1-P^) 
-lOP^ + 2(1-P^) 


2 
2 
2 
2 


Tyr 


CUA 
UUA 
GUA 
AUA 


-IQP^ + 6(1-P^) 

-GP,;^ + i8(i-P(;^) 

-10P(^^ + 6(1-P^) 
-6P^ + 18(1-P(J^) 


2 
2 
-6 
-6 


Ser 


ecu 

UCU 
GCU 
ACU 


6P^ + 2(1-P^) 
-10P(^^ + 2(1-P^) 

6P^ + 2(1-P^) 
-10P($^ + 2(1-P^) 


18 
18 
-6 
-6 


Asp 


cue 

UUC 
GUC 
AUG 


2 

-10P($^ + 6(1-P^) 
2 

-10P(^^ + 6(1-P^) 


18 
18 
-6 
-6 


Glu 


cue 
uuc 

GCA 


2 

-IQP^ + 6(1-P^) 
2 


18 
18 
18 



Continued on next page 
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Continued from previous page 



cl.cL 


aiiticodon 
ACA 


coeff. ch 
-lOP^ + 6(1-P^) 


coeff. Cv 
18 


lie 


CAU 
UAU 
GAU 
AAU 


-lOP^ + 6(1-P^) 
-QP^ + 18(1-P(^) 
-lOP^ + Q{l-P^ 
-QP^ + 18(1-P(^) 


6 
6 

-22 
-22 


Met 


CAU 
UAU 
GAU 
AAU 


-lOP^ + 6(1-P^) 
-6P^ +186(l-Pj) 
-lOPj + 6(1-Pj) 
18Pj + 6(1-P^) 


2 
2 
6 
6 


Asn 


CUU 
UUU 
GUU 
AUU 


-IQP^ + 6(1-P^) 
-QP^ + l^{l-P^) 
-lOP^ + 6(1-P^) 
-6P^ + 18(1-P<J^) 


18 
18 
-6 
-6 


Lys 


CUU 
UUU 
GUU 
AUU 


-lOP^ + 6(1-P^) 
-6Pj + 18(l-Pj) 
-10P<5^ + 6(1-Pj) 
-6Pj + 18(l-Pj) 


18 
18 
18 
18 



Table 6: Value of the coefficients multiplying ch and cy 
in Tav , computed for any anticodon and averaged over 
the two codons for each doublet. 
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